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interactions 



(N 

o 

(N 



Valentin Stanev 1 

1 Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439 

(Dated: May 15, 2012) 

I consider a simple model of a three-band superconductor with repulsive interband interactions. 
The frustration, associated with the odd number of bands, leads to the possible existence of an 
intrinsically complex time-reversal symmetry breaking (TRSB) order parameter. In such state the 
fluctuations of the different gaps are strongly coupled, and this leads to the development of novel 
collective modes, which mix phase and amplitude oscillations. I study these fluctuations using a 
simple microscopic model and derive the dispersion for two physically distinct modes, which are 
gapped by energy less than 2A, and apparently present for all values of interband couplings. 
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I. INTRODUCTION 

Study of the collective modes in superfiuids and super- 
conductors has a long history. In the case of single-band 
systems there are two well-known modes - Bogoliubov- 
Anderson-^ and Schmid 3 , representing oscillations of the 
phase or the amplitude of the order parameter, respec- 
tively. In superfiuids Bogoliubov- Anderson is a gapless 
Goldstone mode, but in superconductors it couples to 
the electromagnetic field and is gapped by energy of the 
order of the plasma frequency. The Schmid mode is al- 
ways gapped by at least 2A (and thus is always damped 
by interaction with quasiparticles). In superconductors 
yet another mode, which couples phase oscillations with 
charge currents, appears close to T c - (it is known as 
Carlson-Goldman). 

In multiband and multicomponent superconductors 
the situation is even richer. This was first realized by 
Leggett j£ who considered a two-band superconductor, in 
which the bands are coupled via a Josephson-like term. 
In such system, apart from the Bogoliubov- Anderson and 
Schmid modes, another collective mode, representing os- 
cillations of the relative phase of the two gaps, is possible 
(observation of this mode in MgB2 has been reported^). 
The Leggett mode is gapped in both superconductors 
and superfiuids (for discussion of Leggett and Carson- 
Goldman modes in a dirty two-band superconductor see 
Ref. 0). There are analogs of this mode in systems 
with two-component order parameters (see, for example, 
Refs. HI andQ), and it has been argued that there is a 
related mode in strongly correlated superconductors 10 . 
This mode exists if the different bands are weakly cou- 
pled and superconducting even in the case in which the 
interband coupling vanishes. In such case there are two 
non-trivial mean-field (MF) superconducting states (with 
phase difference between the two gaps equal to or tt). 
One of these states is metastable (which one depends on 
the sign of the interband coupling). When the interband 
interactions become dominant the metastable state and 
the Leggett mode disappear^ 

Recently, multiband superconductivity has attracted a 
lot of attention in connection with the iron-based high- 
temperature superconductors, since these materials may 
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FIG. 1. The possible superconducting states in a three-band 
system. The angles between the arrows represent the relative 
phase between the different gaps. On the left is the effective 
two-gap s± state, in the middle the three-gap s± state, and on 
the right the TRSB three-gap state (see Ref. [l4| for details). 
The complex state is possible if one or three of the interband 
pairing terms are repulsive. 



have up to four or five superconducting gaps on dis- 
connected pieces of the Fermi surface^. One interest- 
ing feature of these compounds is the likely existence of 
strong repulsive interband interactions. In the case of 
two-band superconductors such interactions lead to an 
unconventional state with a relative minus sign between 
the gaps (the so-called s± state)^. The three-band case 
is even more interesting, since in such system there is an 
intrinsic frustration, and several possible superconduct- 
ing states competed— (schematically shown in Fig[IJ. 
First, the system can effectively behave as a two-gap s± 
superconductor, with one of the bands remaining nor- 
mal. There is also a three-gap s± state, with a n phase 
difference between two of the gaps and the third one. 
The most interesting possibility appears if the three in- 
terband interactions are roughly comparable in strength 
(and thus the two- and three-gap s± states are close in 
energy) - then the system can develop an intrinsically 
complex order parameter, with non-trivial relative phases 
between the different gaps. Such a complex supercon- 
ducting state breaks the time-reversal symmetry. Even 
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if such order parameter is absent in bulk iron-based su- 
perconductors (which are more likely in some variation 
of the s± stated) , a related order parameter can be in- 
duced by bringing in contact iron-based and conventional 
s— wave superconductor a 13 ' 14 . Frustration between the 
s-wave and the s± state at the interface can induce on 
the pnictide side a relative phase, different from the bulk 
7r value. Thus the TRSB state can be stabilized close 
to the boundary, at least under some circumstances (see 
also Ref. M 
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FIG. 2. The two possible phase-amplitude modes around the 
TRSB state. The red/blue/purple arrows represent phase 
oscillations coupled with increase/decrease/no change of the 
amplitude of a given gap. 

In such intrinsically complex TRSB states there is 
a new phenomenon - mixed phase-amplitude collective 
modes. They replace the more conventional Leggett 
modes and appear because phase and amplitude fluctu- 
ations of the gaps on the different bands are coupled 
(unlike the case of multiband s or s± states), due to the 
non-trivial relative phase angle. Recently these modes 
have been discussed in Refs. and I20I using phenomeno- 
logical Ginsburg-Landau (GL) theory™ Within this ap- 
proach it was shown that such modes can have arbitrary 
low mass. In this paper I study the intergap coupling 
using alternative - microscopic - description. This ap- 
proach is very similar in spirit to the one used in Refs. 
Il5l and [l^ . There, however, the coupling between the 
phase and amplitude modes has been mostly neglected. 
In this work I show that such coupling is unavoidable and 
it creates true collective modes, which are gapped by en- 
ergy less than 2 A, and seem to exist for all interband 
interactions. 

Note that coupling between amplitude and phase fluc- 
tuations is mandated by the Galilean invariance even in 
the single-gap case. However, these mixing terms can 
be neglected in the long- wavelength limit.— In contrast, 
the coupling between phase and amplitude fluctuations 
of the different gaps in the TRSB state has a different 
physical origin and survives even in the long- wavelength 



limit. 



II. 



PAIR SUSCEPTIBILITY AND LINEAR 
RESPONSE 



To study the problem I start from a simple generaliza- 
tion of the Bardeen-Cooper-Schrieffer (BCS) model for a 
three-band case: 

ft = J2^- aC U C P a + Yl V P* Cpc^P'a^eq'^ (1) 
p pp'q,a/3 

where a is a band index, the spin indexes are suppressed 
for the moment, and Vp£ is the pairing interaction ma- 
trix, which has both intraband (a = fi) and interband 
(a 7^ 0) components. The TRSB state exists when one 
or three of the interband interactions are repulsive, and 
for concreteness I concentrate on the latter case. The 
intraband interactions can be both negative (attractive) 
or positive (repulsive) , but even then the system can be- 
come superconducting, if the interband interactions are 
sufficiently strong. 

I decouple the interaction terms in a standard BCS 
fashion. This leads to multiband weak-coupling Hamil- 
tonian 

ft = J2^p aC U c P a + a p.« c p +p'« c -p' Q + h - c -'( 2 ) 

P PP',a 

where from the self-consistency condition 

q/3 

(restoring the spin indexes). I will assume that the bands 
are identical - a rather artificial case, but it simplifies the 
calculations considerably (for several examples of non- 
iron based compounds for which this model may apply 
see Ref. M). Thus £ p ,i = £ P)2 = £ p , 3 = £, V 12 = V 13 = 
V 23 = J, V 11 = V 22 = V 33 = V. It is useful to expressed 
A's through the auxiliary variables d: 



dp, a — (Cp+q,a,t c — q,ct ,l) ? 



for which it is straightforward to get: 



di 



J(A 2 



(J + V)A 1 



(J -V)(2J + V) 



(3) 



and analogous expressions for di and d%. 

To study the fluctuations around the TRSB state I 
use standard linear response theory^. In the case of 
identical bands the free energy below T c is minimized by 
the (uniform) state A^mf = A , A 2j a/f = A e 27ri / 3 , 
A3, mf — A e 47 ™/ 3 (up to an arbitrary overall phase)^ 4 -. 
The other mean field solutions - the two and three gap 
s± states - are degenerate and separated from the TRSB 
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state by finite energy. I introduce small fluctuations Here p(r,t) and 9(r, t) represent the amplitude and 
around the TRSB state, and get phase fluctuations, respectively (note that both have di- 

mension of energy). With this the Hamiltonian can be 
Ai = A + pi +i6%, ^ ^ split into two parts -H = Ho+H', where Ho contains 

27ri/3 1 ~ \/3 + i n all the one-particle and mean-field terms, and H! can be 

A2 = A e p 2 



2 2 2l regarded as a (time-dependent) perturbation: 

^oe P3 



n ' = J2 [*p+ q ,i(/'i' 7 i + ^2 + e^7 3 )* P ,i + * P+q 

I <7l+ %/3(T 2 , V3(Tl - fT 2 



-<Tl + \/3<7 2 V^CTi + CT 2 „ , V T 
O P2 5 6*2 + e0<7 3 *p,2 



pq 



*P+q.3 



"P3 



6*3 + e0cr 3 *p. 3 



(4) 



where I used Nambu 2-spinors Vt'p^ = [c p Q ^, Cp^.J, the 
Pauli matrices <7j, momentum- frequency representation 
of 9 and p, and added scalar potential <f>. 

The response of the superconducting state to a pertur- 
bation h Q = [p a , 9 a , ecj>] around its mean field value can 
be written as: 



9 C 



(p > w) = 53nj? a (p,w)^(p ] 



(5) 



where g\ = Tr £ (tft Q a J * p+q , tt ) e - l " t . Using Matsub- 
ara frequencies formalism the polarization tensor is: 

n«,(p,w) = Tr T £ ^~,(i/)(7^~ Pfq ,_ p (w - 1/) 
q,q'^ 

where £7 is the electronic Green's functions in a super- 
conductor: 



toM = -<*£, a *p<,„> = - 



W 2 + |A Q | 2 + £ 2 Q 1 



J p,p 



In writing the linear response in Eq. ([5]) I have used 
the fact that Q is diagonal in band indexes and there 
are no 11^ terms with a ^= (3. This means that the 
coupling between fluctuations in different bands enters 
only through the self-consistency equations. Further- 
more, in the absence of supercurrents, the matrix ele- 
ments coupling amplitude to phase fluctuations of the 
same gap, and to the scalar potential field vanish in the 
long- wavelength limit (Hg a = II 2 ^ = n« = IT^ = for 



each band). Finally, in the model with identical bands 

n«=n« =n«=n« 

With this Eq. ([5]) can be written in a relatively simple 
form. Using the d variables (note that Eq. Q holds for 
each (p, ui) component) the equation for the purely real 
amplitude fluctuations of the first gap can be written as 

2Re(rf 1 (p,w))~n 11 pi(p,w) = 0. 

It is convenient to define 

J „ V 



A 



B 



vo(J- V)(2J + Vy " vo(J-V)(2J + Vy 
and write the expression for d%: 

Re(di) = -u (A + B)Re(Ai) + v Q A(Rc(A 2 ) + Rc(A 3 )) 
= - l y (A + B)p 1 -^-[(p 2 + p 3 ) + V3(9 2 - 9 3 )}, 

where denotes the density of states. 

The equation for the phase fluctuations of Ai is 

2Im(di (p, w)) - n 22 (p, u)0x (p, w) - n 23 (p, u)e4>(p, u) = 0, 

with: 

Im(rfi) = -vq[A + B)9 X + ^[V3( P 2 - p 3 ) - {9 2 + 9 3 )]. 

1 also define Q 23 = U 23 /iy , IP/^o + 

2 J dtl tan( + Ag/2T) / V£ 2 + A 2 = Q u , and ex- 
press the integral through A and B using the gap 
equations (for details see Refs. [27] and |28| ). Utilizing 
these definitions and the equivalence of the bands Eq. 
([5]) can be written in compact matrix form: 
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Q 11 - 2A 


A 


A 





V3A 


-V3A 





A 


Q 11 - 2A 


A 


-y/3A 





V3A 





A 


A 


Q n -2A 


V3A 


-V3A 











-V3A 


V3A 


Q 22 - 2A 


A 


A 


g 23 


V3A 





-V3A 


A 


Q 22 - 2A 


A 


Q 23 


-V3A 


V3A 





A 


A 


Q 22 - 2A 


Q 23 











Q 32 


Q 32 


Q 32 Q 33 - 


-P 2 / 



1 






On 




P3 




9i 












. e 4>. 



= 



The Poisson equation has been used in the last element of 
the lowest matrix row. Note that the coupling constants 
J and V appear in the matrix only as a combination 
in A, which controls all intergap fluctuation couplings. 
This is due to the high symmetries of the model with 
identical bands we are considering. In the case of non- 
identical bands the couplings between the phase-phase 
and phase-amplitude fluctuations are generally different 
and the matrix becomes 19 x 19. 



III. COLLECTIVE MODES 

The collective modes of the system can be obtained 
after analytically continuing to real frequencies. They 
are given by the eigenvectors of the matrix with eigen- 
values equal to zero. The explicit expressions are quite 
unwieldy, but analyzing them reveals a relatively sim- 
ple picture. There are six eigenvectors, and two of them 
represent locked-in oscillations of the phase or amplitude 
of the three gaps, which are trivial generalization of the 
single-gap Bogoliubov- Anderson and Schmid modes. As 
in the single-gap case, the Bogoliubov-Anderson mode 
couples to the electromagnetic field (which gives it a 
mass), while the Schmid mode does not (but is nonethe- 
less gapped by 2Aq). The remaining eigenvectors repre- 
sent the new coupled amplitude and phase fluctuations 
modes (see Fig. There are four eigenvectors and 

two eigenvalues, so each eigenvalue is two-fold degener- 
ate. The energy spectrum of the modes can be obtained 
by solving the equation 

-6A+ Q 11 + Q 22 ± ^/36A 2 + (Q 11 - Q 22 ) 2 = 0. (7) 

For simplicity I consider only the case T = (no ther- 
mally excited quasiparticles). If Q n are small they can 
be expanded to the lowest order in p and i u 27 ' 28 : 

„ ^ 2 ~ °V ~ 4A n 22__^_fV n 23_, 

W ~ 4A 2 ' H ~ 4A 2 ' W 

where c = wf/"\/3 gives the mode velocity. Using these 
results Eq.© simplifies to 

uj 2 = 2A 2 (-6A + 1 ± V36A 2 + l) + c 2 p 2 (8) 

Obviously, the first term on the right side gives the energy 
gap wq of the mixed modes. 




A 



FIG. 3. The energy gap of the mixed modes as function of A, 
as derived from Eqs.© (blue curve) and Eqs.|8} (red dashed 
curve). Apparently the modes can exist both for small and 
large J (negative and positive ^4). 



Before considering the eigenmodes let us look at the 
possible values of A. As seen from its definition, it can be 
positive or negative. When V is finite and attractive (i.e. 
negative) A has infinite discontinuity at J = \V\/2 and is 
negative/positive for smaller/larger values of J. For the 
case in which both interactions are repulsive (positive) 
superconductivity can exist only for J > V, and A is 
always positive and decreasing function of J. 

To obtain ujq these equations have to be solved. Equa- 
tion (UJ) can only be solved numerically, while Eq.© al- 
lows analytic treatment. In both cases it is easy to see 
that the mixed modes are well defined (uq is real and be- 
low 2A ) for positive (negative) A for the plus (minus) 
sign in the equation. Thus for each A there are two dis- 
tinct collective modes (schematically shown on Fig. [5]) 
with identical dispersion. The exact and the approxi- 
mate solutions for ujq are plotted on Fig. [3] and from it 
we see that these modes appear to exist for any given A 
in the interval (—00,00). As the coupling between the 
different gaps decreases (J — > 0, small negative A) ojq 
also becomes smaller. 

Now let me compare these modes with the well-known 
collective modes in one- and two-band superconductors. 
Note that despite the fact that the mixed modes involve 
amplitude oscillations, they are gapped by less than 2A . 
This is quite surprising and is entirely due to the phase- 
amplitude coupling. In contrast, purely amplitude modes 
cannot have energy less then twice the gap. The disper- 
sion of the mixed modes is similar to that of the Leggett 
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mode, but again there are important differences. As the 
results above indicate, the mixed modes are present for 
all values of J, in sharp contrast to the Leggett mode, 
which only exists for weak interband interactions. There 
is one crucial difference between the two- and the three- 
band superconductors, which may explain this - in the 
former the number of the mean-field solutions reduces 
from two to one when the interband interaction starts to 
dominate. In the three-band case such reduction does 
not occur, and there can be three mean-field order pa- 
rameters, even if the intraband interactions are set to 
zero. However, the reader should be warned that using 
the above results in the region of small and positive A 
may be problematic, for at least two reasons. First, this 
region is formally outside of the regime of validity of the 
weak-coupling calculation (since in this region J is large). 
The second and more physical reason is that when in- 
terband interactions dominate it is unclear whether well 
defined relative phase even exist. 

Several features of the above results are due to the 
peculiarity of the particular case with identical bands. 
First, the electromagnetic potential does not couple to 
the mixed modes. In the general case such coupling 
should be expected, since phase fluctuations are conju- 
gate to density fluctuations, and changes in the relative 
density produce charge imbalance, which couples to the 
field. Second, there are two degenerate in energy, but 
physically different modes (shown in Fig. [2]). The de- 
generacy is due to the fact the two metastable s± states 
(states 1 and 2 in Fig. [T]) have the same energy. In the 
case of non-identical bands these states split, and there 
are two different modes with distinct dispersions. 



More detailed and sophisticated calculations are neces- 
sary to study the possible existence of such modes in real 
superconductors, such as iron pnictides. These modes 
can be observed experimentally, for example by Raman 
scatterings^ and their presence can be used as a probe 
of the TRSB state. 



IV. CONCLUSIONS 

In conclusion, I have considered a simple microscopic 
model of a three-band superconductor with repulsive in- 
terband interactions. In such a system the frustration 
associated with the odd number of bands can lead to an 
intrinsically complex TRSB state. In this state two dis- 
tinct collective modes develop, which necessarily mix the 
phase and amplitude oscillations of the different gaps. 
These modes appear to be gapped by energy less than 
2A in the cases of both weak and strong interband cou- 
plings. 
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